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, Single-particle excitation spectra of 7-Fe in the paramagnetic state have been investigated 

by means of the first-principles dynamical coherent potential approximation theory which has 
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recently been developed. It is found that the central peak in the density of states consisting 
of the t2g bands is destroyed by electron correlations, and the Mott-Hubbard type correlated 
bands appear. The results indicate that the 7-Fe can behave as correlated electrons at high 
temperatures. 
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\ The 2>d transition metals are wellknown to behave as a typical itinerant electron system 

in which the ground-state properties are explained weh by the band theory.^ The cohesive 

properties such as lattice parameters and bulk moduli, and the ground-state magnetizations 

are in fact quantitatively explained by the density functional theory for band calculations. 

\ On the other hand, excitations in these metals are often not explained by the band theory. 

CN . The magnetic properties of Fe, Co, and Ni, for examples, show at finite temperatures the 

local moment behaviors as explained by the Heisenberg model. ^ Excitation spectra of Ni 

I observed by means of the X-ray photoelectron spectroscopy (XPS) are wellknown to show 

Q^ \ a d-band narrowing and a satellite peak at 6 eV below the Fermi level, ^ which can not be 

O ■ 

obtained by an independent-particle picture. These results indicate that the effects of electron 
^ correlations in transition metals strongly depends on the details of parameters controlling 

?H ' the physical quantities. Small change in the key parameters might cause anomalous behaviors 

such as high-temperature superconductivity which has recently been found in the iron-arsenide 
system.^ In order to clarify the characteristic features of transition metal systems, one has 
to examine their electronic properties on the basis of realistic band theory which takes into 
account electron correlations. 

In this letter, we present our numerical results of single-particle excitation spectra for 
7-Fe in the paramagnetic state, which are obtained by the first-principles dynamical coherent 
potential approximation (CPA) theory, and demonstrate that 7-Fe can be regarded as a 
strongly correlated electron system at high temperatures, though their ground-state properties 
are believed to be well explained by a band theory. 
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The first-principles dynamical CPA theory^ is tlie dynamical CPA^ combined with the first 
principle tight-binding (TB) linear muffin-tin orbital (LMTO) Hamiltonian.^'^ The former is a 
dynamical version of the single-site spin fluctuation theory developed by Cyrot/^ Hubbard/^ 
and Hasegawa^^ since early in the 1970, and has recently been shown-*^"* to be equivalent to 
the dynamical mean field theory (DMPT).-*^^ The theory describes the electronic and magnetic 
properties at finite temperatures efficiently taking into account the dynamical corrections to 
the spin and charge fiuctuations. Note that unlike the early quantum Montc-Carlo (QMC) 
calculations combined with the DMFT^^ the present approach can treat the transverse spin 
fluctuations for arbitrary d electron number. 

We adopt the TB-LMTO Hamiltonian Hq plus the following intraatomic Coulomb inter- 
actions Hi between d electrons. 

Hi = ^ ^0 num^numi + ^ {Ui - ■^J)niim'hilm' - ^ JSilm ■ him' ■ (1) 

i ITT, m>m' m>m' 

Here Uq (Ui) and J are the intra-orbital (inter-orbital) Coulomb interaction and the exchange 
interaction, respectively. n.iim„ is the number operator for electrons with orbital Im and spin 
a on site i. hum (si/m) with I = 2 is the charge (spin) density operator for d electrons on site 
i and orbital m. 

In the dynamical CPA,^ we transform the interacting Hamiltonian Hi into a dynamical 
potential v in the free energy adopting the functional integral method, and expand the free 
energy with respect to sites after having introduced a uniform medium, (i.e. a coherent po- 
tential) J^Luii^n)- Note that L = (1,171), and u;„ denotes the Matsubara frequency. The first 
term in the expansion is the free energy for a uniform medium, The second term is an 

impurity contribution to the free energy. The dynamical CPA neglects the higher-order terms, 
so that the free energy per atom is given by 

J^CPA = J^P] - /3-HnC J d^e-^^eff(^). (2) 

Here (3 is the inverse temperature. C is a normalization constant. ^ denotes the static field 
variable on a site. -Eeff(^) is an effective potential projected onto the static field ^. It consists 
of the static term E'stC^) and the dynamical correction term £^dyn(^)- The latter is given by a 
Gaussian average of the determinant D of the scattering matrix due to dynamical potential 
as follows. 

g-/3Edy„(^) = det[l - (v - vo)g]. (3) 

Here (v)LriaL'n'a' = VLaa'{i^n - i^n')^LL' ((vq) LnaUn' a' = Vlaa' (0)^LL'^nn') iS the dynamical 

(static) potential, and (g)LnaL'n'a' is the Green function in the static approximation. The 
overline denotes the Gaussian average with respect to the field variables. 

In order to treat the dynamical potential analytically, we expand the dynamical part of 
the effective potential with respect to the frequency v of ^^^^'(^'^i/)) and neglect the mode- 
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mode coupling terms. This is called the harmonic approximation, and D ~ 1 + ^^(i^jy — 1). 
Here is a sub-matrix of D in which the dynamical potential v has been replaced by its i' 
component only. 

The effective medium (i.e., the coherent potential) can be determined by the stationary 
condition 6J-cpa/S^ = 0. This yields the CPA equation as 

{GLait i^n)) = FLaiioJn) ■ (4) 

Note that ( ) at the l.h.s. is a classical average taken with respect to the effective po- 
tential Eefli^), FLa{iuJn) = [(«<^n - Hq - 

iLaiLa is the Coherent Green function, 
where (-ffo)iL(TjL'CT is the one-electron TB-LMTO Hamiltonian matrix, and {'^)iLajL'(7 = 
TjL^{iujn)^ij^LL' ■ Furthermore, GLai^i'^'^n) is the impurity Green function given by 



E 



(5) 



1 + 

Here the first term at the r.h.s. is the Green function in the static approximation. The second 
term is the dynamical corrections, and KLaii^n) = 1 ~ FL^{i(jJn)~'^5FL^{iUn) / 5T,La{i^n)- 

Solving the CPA equation (4) self-consistently, we obtain the effective medium TiL„{iujn). 
The density of states (DOS) for the single-particle excitation spectrum is obtained by means 
of a numerical analytic continuation of the coherent potential. In the numerical calculations, 
we expanded the dynamical correction term -Di/(^) in Eqs. (3) and (5) with respect to the 
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Fig. 1. Density of states (DOS) as the single-particle excitations for 7-Fe in the paramagnetic state 
calculated by the dynamical CPA (solid curve). The DOS in the LDA is shown by dashed curve. 
The BIS^'' and XPS^^ data are shown by dotted curves. 
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Fig. 2. Partial d DOS for eg (dashed curve) and t2g (solid curve) electrons. The total d DOS are 
shown by thin solid curves. The partial DOS in the LDA are shown by dotted curve (eg) and 
dot-dashed curve (t2g). 

Coulomb interactions as ^D^ (^) = 1 +D^^^ (^) +D^^^ (^) + D^f ^ (^) +D^^^ (^) + • • • , and calculated 
the r.h.s. up to the second order exactly as has been made in our previous calculations.^ In 
addition to the second-order terms, the third and fourth order terms are taken into account 
approximately in the present calculations by using an asymptotic approximation.^ We have 
obtained the TB-LMTO Hamiltonian at the observed lattice constant 6.928 a.u. at 1440 K 
using the local density approximation (LDA) in the density functional theory. The average 
Coulomb and exchange energy parameters {U and J) for Fe are taken from the LDA+U 
band calculations^^ ; {U,J) = (0.169,0.066). Uq and Ui are obtained from the relations 
Uo = U + 8 J/5 and Ui = U - 2 J/5. 

We calculated the DOS for 7-Fe in the paramagnetic state at high temperatures (T = 
2000K) where the present theory works best. The results are presented in Fig. 1. The fee DOS 
in the LDA is characterized by the main peak near the top of d bands, the central peak around 
uj = —0.1 Ry, and the third peak near the Fermi level. The first two peaks originate in the 
t2g bands and the third one is due to the eg bands. In the dynamical CPA calculations, the 
central peak is destroyed by electron correlations and the DOS shows the two peaks. The two 
peak structure is more explicitly found in the partial d DOS for Cg and t2g orbitals as shown 
in Fig. 2. Both local DOS show a dip at a; = —0.1 Ry. Especially the central peak of the t2g 
bands in the LDA disappears due to electron correlations and the spectral weight moves to 
the lower and higher energy regions {i.e., uj^ — 0.3 Ry and a; ^ 0.1 Ry). These changes in the 
DOS are caused by a strong scattering peak of —Im.'ELai'^ + around u = —0.1 Ry where 
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the central peak of the t2g band in the LDA is located, and also by a change of the energy 
shift ReSio-('^ + i^) in sign; ReT,L„{LO + i6) < for lu < —0.1 Ry and ReSio-('^ + iS) > 
for u; > —0.1 Ry. The two peaks at a; = —0.25 and lu = 0.05 Ry in the DOS are therefore 
considered to be the lower and upper Mott-Hubbard bands as found in the half-filled Hubbard 
model. 

There is no experimental data for the bulk 7-Fe at high temperatures as far as we know. 
The BIS experimental data for 8 fee Fc monolayers on Cu (100) surface at room temperature^^ 
are consistent with our results as shown in Fig. 1. The photoemission data^^ for 5 fee Fe 
monolayers on Cu (100) are also shown in the figure. The peak at w = —0.205 Ry is usually 
interpreted to be due to the emission from the bulk Cu substrate. Another interpretation 
is that both the lower Hubbard bands for Fe and the Cu-substrate bands are superposed. 
Assuming the latter, the present results are consistent with the experimental data. 

We have also performed the DOS calculations for 7-Mn to make sure the formation of 
the Mott-Hubbard bands in the fee transition metals. We adopted the average Coulomb and 
exchange parameters {U,J) = (0.192,0.061), which arc taken from the LDA -|- U calcula- 
tions with use of the Hartree-Fock-Slater potentials''^ and from the atomic calculations,'^ 
respectively. Wc find again the two-peak structure in the DOS as shown in Fig. 3 when the 
dynamical CPA is applied. The energy difference between the lower peak and the upper one is 
larger than that of 7-Fe by about 0.1 Ry and the valley between the peaks becomes deeper by 
a factor of two. The present results are consistent with those calculated by the QMC-I-DMFT 
with use of the Hamiltonian without transverse spin fluctuations.'' 

There is no XPS experimental data for the bulk 7-Mn. But, the photoemission data" for 
the 20 monolayer fee Mn on the CusAu (100) are consistent with the results of calculations 
as shown in Fig. 3. The high-energy peak at = 0.1 Ry may be justified by the BIS data,'^ 
although the data arc taken for a-Mn at room temperature. 

It should be noted that in the case of Ni the central peak of the t2g bands in the LDA 
shifts to uj = —0.15 Ry. Furthermore, the scattering magnitude — ImS£,a-('^ + i^) is weakened 
by a factor of ten or more as compared with the case of 7-Fe, and it is enhanced around 
u = —0.45 Ry instead of a; = —0.15 Ry. Moreover the energy shift ReSLo-('^ + i6) < for 
UJ > —0.45 Ry and ReT,Lcr{^ + i5) > for u < —0.45 Ry. Therefore the valley found in the 
7-Fe at w = —0.1 Ry disappears in the case of Ni. Instead, a small satellite peak appears 
around a; = —0.50 Ry (6 eV satellite^). 

We note that the Mott-Hubbaxd bands have recently been found in the Fe-based compound 
LaFeAsO. There Fe atoms form a two-dimensional square lattice. The DOS'^ calculated by 
the DMFT clearly show a two-peak structure with a deep depression around u> = —0.1 Ry, 
which is essentially the same as the present results of calculations though the LaFeAsO system 
seems to be considerably stronger than the 7-Fe in Coulomb interaction strength. 
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Fig. 3. The DOS for 7-M11 in the paramagnetic state (sohd curve). The DOS in the LDA is shown 
by dashed curve. The XPS^^ and BIS^^ data are shown by dotted curves. 



In summary, we have calculated the DOS of single-particle excitation spectra for 7-Fe in 
the paramagnetic state on the basis of the first-principles dynamical CPA. We found that the 
7-Fe shows the strong electron correlation effects on the the single-particle excitations at high 
temperature region, i.e., the disappearance of the central peak of the t2g bands due to electron 
correlations and a formation of the Mott-Hubbard type bands. The result is quite different 
from the itinerant electron picture of 7-Fe obtained by the ground-state band calculations.^"^ 
Systematic investigations of the bulk fee transition metals at high temperature region with 
use of the photoelectron spectroscopy are highly desired to verify the present results. 
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